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I review results related to the CKM angle (f>\{0). These results include recent measurements of CP- violation from 
the BaBar and Belle experiments in fe — > ccs, b ccd and b sqq processes. 



1 Introduction 

1.1 The B Physics Program 

The B physics program addresses several funda- 
mental questions. Is the irreducible phase in 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix the 
source of all CP-violating phenomena in the B 
system?^ Or is CP-violation, the first manifesta- 
tion of physics beyond the Standard Model? A re- 
lated question is whether there are new CP- violating 
phases from physics beyond the Standard Model.^ 

The unitarity of the CKM matrix implies the ex- 
istence of three measurable phases. In the convention 
favored at KEK and Belle, these are denoted 



parameters S and A are given by 
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while at SLAC and at BaBar these angles are usually 
referred to as (3, a and 7, respectively. 

As first noted by Bigi, Carter and Sanda,^ there 
are large measurable CP-asymmetries in the decays 
of neutral B mesons to CP-eigenstates. In the decay 
chain T(4S') ^ B^B^ ^ /cp/tag, where one of the 
B mesons decays at time tcp to a final state fcp and 
the other decays at time ttag to a final state /tag that 
distinguishes between P" and 5°, the decay rate has 
a time dependence given by"^ 
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where t^o is the B^ lifetime, Ato^ is the mass dif- 
ference between the two P° mass eigenstates. At = 
tcp — ttag, and the 6-flavor charge q = +1 (—1) when 
the tagging B meson is a P" (P")- The CP- violation 
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where A is a complex parameter that depends on 
both the B'^B^ mixing and on the amplitudes for P° 
and P'^ to decay to fcp- To a good approximation, 
the SM predicts S = — ^/ sin2(/)i, where ^/ = 1) 
corresponds to CP-even (-odd) final states. Direct 
CP-violation, ^ = (or equivalently |A| = 1), is 
expected for both h ccs and h sss transitions. 

1.2 Accelerators and Detectors 

The P-factory accelerators, PEPII'' and KEKB^ 
were commissioned with remarkable speed starting 
in late 1998. The experiments, BaBar^ and Belle,"^ 
started physics data taking in 1999. In the summer 
of 2001, the two experiments announced the obser- 
vation of the first statistically significant signals for 
CP-violation outside of the kaon system.*'^ 

Due to the extraordinary performance of the two 
accelerators, the most recent results reported in the 
summer of 2003 at the Lepton-Photon Symposium 
are based on very large data samples. BaBar has in- 
tegrated 113 fb~^ on the T(AS) resonance while Belle 
has integrated a sample of 140 fb~^. KEK-B also 
passed a critical milestone for e^e" storage rings and 
achieved a peak luminosity above 1 x 10"^'* cm~^s~^. 

1.3 The Principle of the Measurement 

The measurement of time-dependent CP-asymmetry 
requires: 

• A large sample of T(4S') decays into B^B^ pairs. 
To boost the T(4S') decay frame so that the 
P mesons' flight length can be measured with 
solid-state vertex detector technology, both the 
KEKB and PEP-II accelerators use asymmetric 
energy beams with energies of 8.0 and 3.5 GeV 
or 9.0 and 3.1 GeV, respectively. 
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Figure 1. The fully reconstructed CP-eigenstate sample used 
by Belle. This sample is obtained from a data sample with an 
integrated luminosity of 140 fb~^. 
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Figure 2. The Pg {B momentum in the CM frame) distribu- 
tion for the B — > J/ipKi^ sample used by Belle. This sample 
is obtained from a data sample with an integrated luminosity 
of 140 fb~^. The shaded portions show the contributions of 
different background components. The vertical dashed lines 
indicate the signal region. 



• EfRcient reconstruction of B ^ X^cK^ de- 
cays. This implies accurate measurements of 
momenta and energies of neutrals using CsI(Tl) 
crystal calorimeters in addition to good charged 
particle tracking in small cell drift chambers and 
efficient identification of leptons and as well 
as mesons. 

• A measurement of At. This is related to the 
measurement of Az, the spatial distance be- 
tween the decay vertices and achieved at both 
experiments by using double-sided silicon strip 
detectors situated at small radii close to the in- 
teraction point. 

• A determination of the flavor of the accompa- 
nying B ("tagging"); this is based on the iden- 
tification of electrons, muons and charged kaons 
and the measurement of their charge. 

More detailed descriptions of the detectors®'^ 
and the experimental analysis procedure are avail- 
able elsewhere.^° 

2 Status of CP- Violation in 6 ^ ccs 
Processes 

Belle and BaBar reconstruct B^ decays to the fol- 
lowing b — > CCS CP-eigenstates: J/4>Ks, iIj{2S)Ks, 



XciKs, ilcKs for ^/ = -1 and JlitKL for ^/ = 
+1." The two classes (^/ = ±1) should have CP- 
asymmetries that are opposite in sign. 

Both experiments also use P" — » J/tpK*'^ decays 
where K*'^ Kst^^ ■ Here the final state is a mixture 
of even and odd CP. The CP content can, however, 
be determined from an angular analysis of other il)K* 
decays. The CP-odd fraction is found to be small 
(i.e. (19 ± 4)% ((16 ± 3.5)%) in the Belle (BaBar) 
analysis) . 

The most recent BaBar analysis is based on a 
data sample with an integrated luminosity of 81 fb~^ 
and was first presented in 2002.^ There is a corre- 
sponding published Belle result also shown in 2002 
with 78 fb^^.^ At this Symposium, Belle provided a 
new preliminary result for their 140 fb^^ sample.^^ 

The data sample used for the recent Belle mea- 
surement is shown in Fig. 1 and Fig. 2. Table 1 
lists the numbers of candidates, A^ev, and the esti- 
mated signal purity for each fcp mode. It is clear 
that the CP-eigenstate samples that are used for the 
CP- violation measurements in & — > ccs are large and 
clean. 

In the summer of 2001, the first statistically sig- 
nificant measurements of the CP- violating parame- 
ter sin 2(/)i were reported by Belle and BaBar. Belle 



Table 1. The yields from Belle for reconstructed B — > fcp 
candidates after flavor tagging and vertex reconstruction, A^ev, 
and the estimated signal purity, p, in the signal region for each 
fcp mode. J/tp mesons are reconstructed in J/xp — > p.'^ fJ.~ 
or e+e" decays. Candidate mesons are reconstructed in 
K'i — > 7r+7r~ decays unless otherwise written explicitly. 
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-1 


1997 


0.976 ± 0.001 


J/iPKliTrOnO) 


-1 


288 


0.82 ±0.02 


■>p{2S){e+e-)K0 


-1 


145 


0.93 ±0.01 


ii{2S)(J/i>TX+-K-)Kl 


-1 


163 


0.88 ±0.01 




-1 


101 


0.92 ±0.01 


r),{KlK--K+)K° 


-1 


123 


0.72 ±0.03 


r)c{K+K--K°)Kl 


-1 


74 


0.70 ±0.04 


ric{pp)K% 


-1 


20 


0.91 ±0.02 


AU with (,j = -1 


-1 


2911 


0.933 ± 0.002 


J/Vii'*"(i^iJ7r") 


+1(81%) 


174 


0.93 ±0.01 




+1 


2332 


0.60 ±0.03 



found 

sin 201 = 0.99 ±0.14 ±0.06 (5) 

while BaBar obtained 

sin 201 = 0.59 ±0.14 ±0.05. (6) 

The resuhs were based on data samples of compara- 
ble size (31 million and 32 million BB pairs, respec- 
tively). 

The new Belle data are shown in Fig. 3. This 
figure shows the At distributions where a clear shift 
between B^ and B^ tags is visible as well as the 
raw asymmetry plots in two bins of the flavor tag- 
ging quality variable r. For low-quality tags (0 < 
r < 0.5), which have a large background dilution, 
only a modest asymmetry is visble while in the sub- 
sample with high quality tags (0.5 < r < 1.0), a very 
clear asymmetry with a sine-like time modulation 
is present. The final results are extracted from an 
unbinned maximum-likelihood fit to the Ai distribu- 
tions that takes into account resolution, mistagging 
and background dilution. The new Belle result with 
140 fb-i (152 million BB pairs) is 

sin 201 = 0.733 ± 0.057 ±0.028. (7) 

The new Belle result may be compared to the 
BaBar result with 78 fb^^ of 

sin 201 = 0.741 ±0.067 ±0.03. (8) 

Both experiments are now in very good agree- 
ment. A new world average can be calculated from 
these results, 

sin 201 = 0.736 ±0.049. (9) 
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Figure 3. Belle data from 2003: (a) Ai distributions for B° 
and B'^ tags (b) raw asymmetry for low-quality tags and (c) 
raw asymmetry for high-quality tags. The smooth curves are 
projections of the unbinned likelihood fit. 

This world average can be interpreted as a constraint 
on the CKM angle 0i. This constraint can be com- 
pared to the indirect determinations on the unitarity 
triangle.^'^ This comparison is shown in Fig. 4 and is 
consistent with the hypothesis that the Kobayashi- 
Maskawa phase is the source of CP- violation. 

The measurement of sin(20i) in 6 ^ ccs modes, 
although still statistically limited, is becoming a pre- 
cision measurement. The systematics are small and 
well-understood. Recently, BaBar physicists discov- 
ered a new small source of systematic uncertainty 
due to CP-violation in h ^ cud decays on the tag- 
ging side.^"* 

The presence of an asymmetry with a cosine 
dependence (|A| ^ 1) would indicate direct CP- 
violation. In order to test for this possibility in 
h CCS modes. Belle also performed a fit with 

acp = ^O-'^^'^/l'^l ^^"-^ l-^l ^^^'^ parameters, keep- 
ing everything else the same. They obtain 

|A| = 1.007 ±0.041(stat) (10) 
acp = 0.733 ±0.057(stat). 
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Figure 4. Indirect constraints on the angles of the CKM tri- 
angle compared to the most recent direct measurements of (/>i 
from Belle and BaBar. The theoretical uncertainties in the 
indirect constraints are conservatively estimated by the CKM 
fitter group. 



for all the b — > ccs CP modes combined. This re- 
sult is consistent with the assumption used in their 
primary analysis. 



Studies of CP- Violation in b 
Processes 



ccd 



Neutral B decays to CP-eigenstates that proceed by 
b — > ccd processes are expected to have the same CP- 
violation as 5 — > CCS since both are sensitive to the 
phase of B — B mixing. A small deviation from this 
expectation is possible because of the contribution of 
b ^ d penguin diagrams (a.k.a. "penguin pollution" ) 
in the decay modes that are examined. Penguin pol- 
lution may also give rise to direct CP-violation and a 
CP-violating term with a cos(Ar7ic;Ai) dependence. 

The b ccd decay modes that have been used 
so far for CP- violation studies are B — + D*^D*^ , 
B D*+D-, and B J/i^TT^}^-^^ The effect of 
penguin pollution might be expected to be the largest 
in B J/ipiT^ because the penguin contribution is 
not color-suppressed in that mode. 

For B ^ ipTr°, with 81 ib'^ BaBar has a signal 
of 40 ± 7 events and finds 

sin20ie//(B ^ 4>n") = 0.05 ± 0.45 ± 0.16. (11) 



The corresponding result from Belle is based on 140 
fb^^ and uses 89 ± 10 events.-'^^ They obtain 

sin20ie//(P ^ V7r°) = 0.72i|]-j^ ± 0.08. (12) 

In both cases, the systematic error includes the pos- 
sibility of CP- violation in a small component of the 
background that peaks under the signal. 

The b — > ccd mode B — > D*^D*~ has a vector- 
vector final state and requires special treatment since 
it includes contributions from both CP-even and odd 
components. To extract the CP-odd fraction, one 
fits the angular distribution in the transversity basis. 
The result from BaBar based on a sample with 156± 
14 signal events is, 



R±_ = 0.063 ± 0.055 ± 0.009, 



(13) 



where the quantity R±^ is the fraction of the CP- 
odd component. The measurement indicates that 
B° D*+D*- is mostly CP-even. 

The time distributions from BaBar for B — > 
£)*+£)*- are shown in Fig. 5. BaBar finds 

sin20ie//(B D*+D*-) = -0.05 ± 0.29 ± 0.10, 

(14) 

which is about 2.5cr from the result in 6 — s- ccs modes. 
This may be a statistical fiuctuation or could be an 
indication that the Standard Model b ~* d penguin 
contribution is large. The fit includes the possibility 
of direct CP- violation. The parameter A is found to 
be 0.75 ± 0.19 ± 0.02, which is consistent with unity, 
as expected for no direct CP-violation. 

Since P*^ D*^D^ and its charge conju- 
gate are not CP-eigenstates, a modified treatment 
is required. There are four rather than two CP- 
violating observables that are determined from a 
time-dependent fit to the different D*D charge 
states. 

BaBar finds, 

S'_^_ = -0.82 ±0.75 ±0.14, (15) 

5_+ = -0.24 ±0.69 ±0.12, (16) 

= +0.47 ±0.40 ±0.12, (17) 

= +0.22 ±0.37 ±0.10. (18) 

In the limit of no penguins and assuming factor- 
ization in these hadronic decays, S y = = 

— sin 2(j)i and = A ^ = 0. The above results 

for CPV in B D*D decays are consistent with 
this limit. 
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Figure 5. BaBar results on CP-violation in _B — > D*+D*~. 
The top two figures show the At distributions for and 
tags. The third plot shows the raw time asymmetry distribu- 
tion. 

Observation of the CP-eigenstate mode B — > 
D^D~ was reported by Belle at this conference. 
With 140 ih-\ the 5cr signal contains 24.3 ± 6.0 
events. In the future, this mode can also be used for 
time-dependent measurements of CPV in & ^ ccd 
processes. 

The results of CP- violation measurements for 
h — > ccd decays are summarized in Fig. 6. The mea- 
surements are not yet precise enough to definitively 
demonstrate the presence of penguin pollution. 

4 Status of CP- Violation in 6 ^ sqq 
Penguin Processes 

In addition to the program of measuring the other 
remaining angles of the unitarity triangle that is dis- 
cussed in the contribution by Jawahery,^^ there is 
also the question of whether there are additional CP- 
violating phases from new interactions or physics be- 
yond the Standard Model. At the moment, such new 
phases are poorly constrained. 

One way to attack this question is to mea- 
sure the time-dependent CP-asymmetry in penguin- 
dominated modes such as P° — > (fiK^, ^ -q 
or B'^ Kgir'^, where heavy new particles may con- 
tribute inside the loop, and compare it to the asym- 
metry in — » J/ipKg and related b ccs charmo- 
nium modes. 



D AD A □ no 

-0.05±0.49±0.16 1 • 

°£ Belle 03 

^ 0.72 ±0.42 ±0.08 

Average 

0.40 ±0.33 


1 

1 

1 — * — 


i 


1 


b BABAR 03 

+ -0.06±0.37±0.13 1 » 

Q 


1 




BABAR 03 [+-] 

1 0.82±0.75±0.14 
Q 

+ BABAR 03 [-+] 


1 


■ 1 


Q 0.24±0.69±0.12 1 




if 




Charmonium (average) 
0.736 ±0.049 










-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 



sin(2P3„) 

Figure 6. Summary plot of results on CP-violation in fe ^ ccd 
modes. 




Figure 7. Beam constrained mass and AE distributions for 
B K°g-K° from BaBar. 

The mode B Ksi^^ proceeds through a 6 ^ 
sdd transition. The BaBar data on P ^ K^tt^ are 
shown in Fig. 7. To be useful for time-dependent 
CPV studies at least one of charged pions from the 
K'g must be detected in the BaBar silicon vertex 
detector.^° There are 123 ± 16 events of this type 
that are then used to obtain 

sin2,/)ie//(P^ifg7rO) =G.48iEJ-j8±0.11. (19) 

The time distributions are shown in Fig. 8. The di- 
rect CP-violation parameter is A = —0.40^0 27 ^ 
0.10.^° When A is fixed to zero, the value of S' = 
sin(2(/>ie//) shifts slightly to 0.411°;^^ ±0.11. The 
results for B — > Kstt^ are consistent with the value 
from the b — > ccs modes, sin2(/)i = 0.736 ± 0.049. 

The mode B v'^s expected to include con- 
tributions from b — > suu and b — > sdd penguin pro- 
cesses. The beam constrained mass distribution for 
the B v'^s sample used by Belle is shown in 



6 




At (ps) 

Figure 8. BaBar data on B K°TT°. The top two fi eures 
show the At distributions for and tags, separately. 
The third plot shows the raw time asymmetry distribution. 



Fig. 9 and contains 244 ± 21 signal events.^^ Belle 
finds (Fig. 10), 

sin20ie//(B ^ rj'K"^) = 0.43 ± 0.27 ± 0.05 (20) 

The BaBar data is shown in Fig. 11. They obtain, 

sin 20ie//(B ^ Tj'K"s) = 0.02 ± 0.34 ± 0.03 (21) 

The average of these two results for B v'^s 
about 2.2a from the b —>■ ccs measurement, which is 
the Standard Model expectation. 





5.22 5.24 5.26 

(QeV/c^) 



5.22 5.24 5.26 

(QeV/c^) 



Figure 9. Beam constrained mass distributions for B v' 
(left) and B K+ K- (right). 



The decay mode B K+R-Rl, where K+R- 
combinations consistent with the (j) have been re- 
moved, is found by Belle to be dominately CP-odd^^ 
and thus can be treated as a CP-eigenstate and 
used for studies of time-dependent CP- violation in 




Figure 10. Belle data for the raw asymmetry in B^ —> rj'Kg. 
The upper plot shows the data for low-quality tags while the 
lower plot shows the higher quality tags. The dashed curves 
are the expectations from the Standard Model. 



b — > sqq processes. The beam constrained mass dis- 
tribution for the B K~^K~Kg sample used by 
Belle is shown in Fig. 9. There are 199 ± 18 signal 
events. Belle obtains, 



-K"g) = 0.51±0.26±0.05tE;j^, 



sin2<?!)ie//(P K+K 

(22) 

where the third error is due to the uncertainty in 
the CP content of this final state.^^ The results for 
B K^K^Kg are also consistent with b ccs 
decays. However, in this decay there is also the pos- 
sibility of "tree-pollution" , the contribution of the 
b — > uus tree amplitude that may complicate the 
interpretation of the results.^^ 

The P° decay, which is dominated by 

the b — > sss transition, is an especially unambiguous 
and sensitive probe of new CP- violating phases from 
physics beyond the SM.^'' The SM predicts that mea- 
surements of CP-violation in this mode should yield 
sin2(/)i to a very good approximation.^^'^'^ A signifi- 
cant deviation in the time-dependent CP-asymmetry 
in this mode from what is observed in 6 — > ccs decays 
would be evidence for a new CP- violating phase. 

The B — > 4>Kg sample used by BaBar is shown 
in Fig. 12. The signal, obtained from a sample with 
an integrated luminosity of 110 fb~^ , contains 70 ± 9 
events.^" The time distributions for the BaBar data 
are shown in Fig. 13. They obtain 

sin2</)ie//(P ^ (t)K°s) = 0.45 ±0.43 ±0.07. (23) 
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Figure 11. BaBar data on i? — > r/'K^. The top two figures 
show the At distributions for and tags, separately. The 
third plot shows the raw time asymmetry distribution. 




Figure 13. BaBar time difl'erence and asymmetry data dis- 
tributions in B — > (pK^. The top two figures show the At 
distributions for and B^ tags, separately. The third plot 
shows the raw time asymmetry distribution. 




Figure 12. Beam constrained mass distributions for B — > if>K\ 
from BaBar (left) and Belle(right). 
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Figure 14. Belle data: (left) At distributions for low-quality 
tags and (right) for high-quality tags. The dashed curves show 
the background contributions. 



This value is consistent with the Standard Model 
expectation, but is somewhat different from the 
value obtained with the 81 fb~^ sample, which was 
sin2(/)ie// = -0.18 ±0.51 ±0.09. The new result in- 
cludes more data and a reprocessing of the old data 
sample. After extensive checks with data and Toy 
Monte Carlo studies, the large change in the central 
value is attributed to a la statistical fluctuation.^^ 

The B — > (f>Kg sample used by Belle is shown in 
the right panel of Fig. 12. The selection criteria are 
described in detail elsewhere.^^'^^ The signal contains 
68 ± 11 events. Figure 15 shows the raw asymme- 
tries from Belle in two regions of the flavor-tagging 
parameter r. While the numbers of events in the two 
regions are similar, the effective tagging efficiency is 



much larger and the background dilution is smaller 
in the region 0.5 < r < 1.0. The solid curves show 
the results of the unbinned maximum-likelihood fit 
to the At distribution. 

The observed CP-asymmetry for — ^ 'pKg in 
the region 0.5 < r < 1.0 (Fig. 15 (lower panel)) indi- 
cates the difference from the SM expectation (dashed 
curve). Note that these projections onto the At axis 
do not take into account event-by-event information 
(such as the signal fraction, the wrong tag fraction 
and the vertex resolution) that is used in the un- 
binned maximum likelihood fit. 

The contamination of K^K^Kg events in the 
(j>K'^ sample (7.2 ± 1.7%) is small. Finally, back- 
grounds from the — > fo{98Q)Kg decay, which has 
the opposite CP-eigenvalue to (j)Kg, are found to be 
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Figure 15. Belle data for the raw asymmetry in — > (ji^s- 
The upper plot shows the data for low-quality tags while the 
lower plot shows the higher quality tags. The dashed line is 
the expectation from the Standard Model. 



small (1.6l}'g%). The influence of these backgrounds 
is treated as a source of systematic uncertainty. 
Belle obtains 



sin20ie//(B ^ (l)K[ 



-0.96 ±0.5 



4-0.09 
-0.11 



(24) 



from their likelihood fit to the 4>Kg data. The like- 
lihood function is parabolic and well-behaved. An 
evaluation of the significance of the result using the 
Feldman- Cousins method and allowing for system- 
atic uncertainties shows that this result deviates by 
3.5cr from the Standard Model expectation.^^ 

The Belle group performed a number of val- 
idation checks for their B — > <t>Kg CP-violation 
result. Fits to the same samples with the di- 
rect CP- violation parameter A fixed at zero yield 
sin2(f>ieff = -0.99 ± 0.50(stat) for B° (j)Kl. As 
a consistency check for the S term, the same fit pro- 
cedure is applied to the charged B meson decays 
B+ (l)K+. The resuh is 5 = -0.09 ± 0.26(stat), 
A = +0.18 ± 0.20(stat) for B+ (j)K+ decay. The 
results for the S term is consistent with no CP- 
asymmetry, as expected. The asymmetry distribu- 
tion is shown in Fig. 16. In addition, the (j)K^ side- 
band has been examined as shown in Fig. 16. No 
asymmetry is found in that sample. 
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Figure 16. Belle data: consistency checks of the B 'P-^s 
analysis. The asymmetries in (a) the B^ <j>K^ sample and 
(b) the B — > <f>K'^ sideband sample. 
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Figure 17. Summary plot of results on sin2(/>i and sin2</<ie// 
in fe — ► CCS and b sqq modes. 



5 Conclusion 

Belle presented a new measurement of time- 
dependent CP-violation in b ccs CP-eigenstates. 
This result and previous results from BaBar are in 
good agreement with each other and with the hy- 
pothesis that the Kobayashi-Maskawa phase is the 
source of CP-violation. 

Studies of CP-violation in 5 — > ccd modes are 
progressing. In B D*^D*^ decays, BaBar ob- 
serves a 2.5(7 hint for penguin pollution. More data 
and measurements are needed to clarify whether pen- 
guin pollution is present in this class of decays. 

In P ^ 4>K'g decays there was a surprise. With 
140 fb^^ Belle observed a 3.5(j deviation from the 
Standard Model expectation. This could be an in- 
dication of new physics from heavy particles in the 
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b ^ ss penguin loop. However, BaBar's value moved 
closer to the Standard Model with the addition of 
new data and reprocessing. More precise measure- 
ments of the other b sqq modes can further con- 
strain phases from new physics. For example, new 
physics may contribute differently to pseudoscalar- 
vector and pseudoscalar-pseudoscalar modcs.^ 

The results of CP-violation measurements for 
b sqq penguin decays are summarized in Fig. 17. 
The world average for all 6 — > s penguin decays 
(shown by the dotted line) appears to be displaced 
from the average for b — » ccs modes. The high en- 
ergy physics community will require that this exper- 
imental issue be resolved conclusively in the future. 
This will require large data samples with integrated 
luminosities of at least 1 ab~^ or 1000 fb~^. 
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DISCUSSION 

Stefan Spanier (University of Tennessee): 

1) Unfortunately, the plenary session gives the 

audience only a limited chance to help you to es- 
tablish the results by asking detailed questions. 

2) Knowing the previous value of 5* = —0.7 ±0.6 
from Belle, the newly added statistics must lead 
to an unphysical value of S' < —1.4 leading typi- 
cally to large correlations in S and C (pathologi- 
cal behavior) in this new sample. How probable 
is the value in the new sample? 

3) How strong is the CP-asymmetry in the back- 
ground? 

Tom Browder: 

1) A special breakout session is planned later in 

the Symposium. 

2) For a true value near S = —1, the values 
in the new sample are quite consistent with 
Toy Monte Carlo studies. There is no statisti- 
cally pathological behaviour in either old or new 
data samples. The observed errors are actually 
slightly larger than expected. 



3) The background from B — > foKg and 
B K+K-Kl decays is small and the CP- 
asymmetry from these backgrounds is included 

in the systematic error. 

Alex Kagan (Cincinnati): Can you show the raw 
BaBar data for S{(f)K^) again? 

Tom Browder: Yes. Note that a figure with this 
data was included in the talk and appears in 
the Proceedings. 

Hitoshi Murayama (Berkeley): On the <^-f^5 

mode, the change in the BaBar result was at- 
tributed to a statistical fluctuation. They have 
added only 40% more data. How is that possi- 
ble? Do you have a breakdown of the asymme- 
try between the previous and new data samples? 

Tom Browder: Not only was more data added, but 
the old BaBar data sample was also reprocessed. 
After reprocessing, a small number of events 
changed from P° tags to P° tags (or vice versa). 
This accounts for the shift in the central value. 
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We present the first study of the time-dependent CP- violating asymmetry in B° J/ip tt" decays 
using e'^ e~ annihilation data collected with the BABAR detector at the T(4S) resonance during 
the years 1999-2002 at the PEP-11 asymmetric-energy B Factory at SLAC. Using approximately 88 
million BB pairs, our results for the coefficients of the cosine and sine terms of the CP asymmetry 
are Cj^^^o = 0.38 ± 0.41 (stat) ± 0.09 (syst) and Sj^^^o = 0.05 ± 0.49 (stat) ± 0.16 (syst). 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 

The Standard Model of electroweak interactions de- plex phase in the three-generation Cabibbo-Kobayashi- 
scribes CP violation in 5-meson decays by a com- Maskawa (CKM) quark-mixing matrix [1]. The h^ccs 
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modes such as B° J/ip K'^ yield precise measurements 
of the quantity sin2/3, where /3 = Siig[—V^^V*f^/Vf.^Vi-l] 
(see for example Refs. [2-4]). The decay B° J/ipn^ is 
a Cabibbo-suppressed h ccd transition. In the Stan- 
dard Model both B° J/ip K° and B° J/ip Tr° have 
penguin amplitudes with the same weak phase as the 
tree amplitude, and an additional penguin amplitude 
with a different phase. In B° J/ipK^, the penguin 
amplitude with a different weak phase is suppressed by 
^CKM' where Xckm is the sine of the Cabibbo angle, 
while in B° J/ipn^, the tree and each penguin am- 
plitude are equal to leading order in Xckm- Therefore, 
B° J/ip Tr° may have a CP asymmetry that differs 
from that of B° J/ip K^, with the size of the asymme- 
try serving as a probe of the penguin decay amplitudes 
in both modes. 

BABAR has previously measured the B° J/ip 
branching fraction, (2.0±0.6 (stat)±0.2 (syst))xlO-5 [5], 
using T{4:S) BB decays. For the CP asymmetry mea- 
surement, the flavor (5*^ or 5*^) of the B meson that 
decays to J/ipir'''' is inferred, or tagged, using properties 
of the other B meson and the time evolution of the BB 
system. The decay time distributions, f_|_(f_), of B de- 
cays to a CP eigenstate with a 5*^ (-B*^) flavor tag, are 
given by 



f±(At) 



-lAtl/r^o 



l±5j/^^o sin(AmdAt) 



TCj/^tt" cos {AnidAt) 



where At = tree— ^tag is the difference between the proper 
decay time of the reconstructed B meson and the proper 
decay time of the tagging B meson, Tgo is the B° life- 
time, and Anid is the B°-B° oscillation frequency. The 
coefficients can be expressed in terms of a complex pa- 
rameter A, which depends on both the B°-B° oscilla- 
tion amplitude and the B° and B° decay amplitudes 
to this final state [6]: Sj/^j^o = 2XmX/{l + |Ap) and 
Cj/^^o = (1 — |Ap)/(l -|- |Ap). A decay amplitude with 
only a tree component would give Sjj^T^o = — sin2/3 and 



C 



0. 



The data used in this measurement were collected with 
the BABAR detector [7] at the PEP-II storage ring in the 
years 1999 to 2002. Approximately 81 fb"' of e+e" anni- 
hilation data recorded at the T(45') resonance are used, 
corresponding to a sample of approximately 88 million 
BB pairs. An additional 5fb~^ of data collected approx- 
imately 40 MeV below the T(45') resonance are used to 
characterize non-55 background sources. 

B° J/ip Tr° candidates are selected (details are given 
in Ref. [5]) by identifying J/ip e+e" or J/ip IJ-^ IJ-~ 
decays and tt*^ 77 decays. For the J/ip e+e" {J/i' 
— 7- jj,'^ IJ,~) channel, each lepton candidate must be consis- 
tent with the electron (muon) hypothesis. The invariant 



TABLE I: Efficiencies for the requirement on the Eisher dis- 
criminant and flavor tagging, given independently, with sta- 
tistical uncertainties. 



Type of event 


Efliciency (%) 




Eisher 


Tagging 


B° ^ J/V- tt" 


99.2 ± 0.1 


65.6 ± 0.6 


B° J/4)K°s{n°n°) bkg. 


98.9 ± 0.1 


65.6 ± 0.6 


Inclusive J/tp bkg. 


94.9 ± 0.7 


70.4 ± 1.4 


BB generic bkg. 


98.5 ± 0.4 


61.1 ± 1.6 


Continuum bkg. 


28.6 ± 0.7 


52.3 ± 0.8 



mass of the lepton pair is required to be between 2.95 and 
3.14GeV/c2, and 3.06 and 3.14GeV/c2, for the electron 
and muon channels, respectively. The photon candidates 
used to reconstruct the tt° candidate are identified as 
clusters in the electromagnetic calorimeter (EMC) with 
polar angles between 0.410 and 2.409 rad, that are spa- 
tially separated from every charged track, and have a 
minimum energy of 30 MeV. The lateral energy distribu- 
tion in the cluster is required to be consistent with that 
of a photon. The invariant mass of the photon pair is re- 
quired to between 100 and 160MeV/c^. Finally, the J/ip 
and Tr° candidates are assigned their nominal masses and 
combined using 4-momentum addition. 

Two kinematic consistency requirements are applied 
to each 5 candidate. The difference, A5, between the 
5-candidate energy and the beam energy in the e+e" 
center-of-mass (CM) frame must be —0.4 < A5 < 
0.4 GeV. The beam-energy-substituted mass, mEs = 
^(^/2)2 — (p^)2, must be greater than 5.2GeV/c^, 
where ^/s is the total CM energy and pg is the 5- 
candidate momentum in the CM frame. 

A linear combination of several kinematic and topo- 
logical variables, determined with a Fisher discriminant, 
provides additional separation between signal and e+e" 
qq (q = u, d, s, c) continuum background events. The 
Fisher discriminant uses the following inputs: the zeroth- 
and second-order Legendre polynomial momentum mo- 
ments {L, = Ip*! and = Ip*! where 
p* are the CM momenta for the tracks and neutral 
calorimeter clusters that are not associated with the sig- 
nal candidate, and 6i are the angles between p* and 
the thrust axis of the signal candidate); the ratio of 
the second-order to zeroth-order Fox- Wolfram moments, 
again using just tracks and clusters not associated with 
the signal candidate; |cos6't|, where 6t is the angle be- 
tween the thrust axis of the 5 candidate and the thrust 
axis of the remaining tracks and clusters in the event; 
and I cos 61 \ , where 61 is defined as the angle between the 
negative lepton and 5 candidate directions in the J/ip 
rest frame. The requirement placed on the Fisher dis- 
criminant is 99% efficient for signal and rejects 71% of 
the continuum background. The efficiencies for satisfy- 
ing this requirement are summarized in Table I. 
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We split the backgrounds into four mutually exclu- 
sive categories, two of which have a J/^ from B de- 
cays (B J/ip X). The first background category is 
B° J/ip K^^TT^TT^) decays where one of the mesons 
is nearly at rest in the e+e" CM frame. The second 
background category consists of other B J/ip X de- 
cays (inclusive J/ip ) , which contribute through random 
combinations of J/ip and candidates. The third and 
fourth categories consist of random combinations of par- 
ticles in BB decays (B B generic) and continuum events, 
respectively. Monte Carlo simulation [8] is used to model 
aspects of the B° J/ip li'^^n^TT^), inclusive J/ip , and 
BB generic backgrounds. A sample (■//V'fake) selected 
from data taken below the T(45') resonance is used to 
model the continuum background. In this case, the J/ip 
candidate is reconstructed from two tracks that are not 
consistent with a lepton hypothesis. Monte Carlo simula- 
tion is used to check that this procedure, which increases 
the size of the sample, correctly models the continuum 
background. 

The algorithm for 5-fiavor tagging assigns events to 
one of four hierarchical, mutually exclusive tagging cat- 
egories, and is described in detail in Ref. [3]. The to- 
tal tagging efficiency for the signal and each background 
source is given in Table I. Untagged events are excluded 
from further consideration. Vertex reconstruction and 
the determination of At follow the techniques detailed in 
Ref. [9]. We require —20 < At < 20 ps and an estimated 
uncertainty on At of less than 2.4 ps. 

We extract the CP asymmetry by performing an un- 
binned extended maximum likelihood fit. The likeli- 
hood is constructed from the probability density func- 
tions (PDFs) for the variables mEs, AE, and At. The 
quantity that is maximized is the logarithm of 



N\ 



N 5 

HE 

i=i j=i 



/f n -pf 



(2) 



where 



the number of events for each of the five 



hypotheses (one signal and four background) and N is 
the number of input events. The Vj' are the one- or two- 
dimensional PDFs for variables d, for each signal or back- 
ground type. The parameters are the tagging frac- 
tions for each of the tagging categories a,- (assigned for 
each event i) and each of the signal or background types 
j. For the B° J/iP Tr° signal and B° J/iP K°{Tr°Tr°) 
background, the values of are measured with a sam- 
ple (5flav) of neutral B decays to fiavor eigenstates con- 
sisting of the channels D''*^~ [h~^ = tt~^ , p'^ , and a^) 
and J/ipK*°{K*° A'+tt") [3]. Monte Carlo simulation 
is used to estimate the values for the inclusive J/ip 
and BB generic backgrounds, while the J/'tpfg^^ sample 
is used for the continuum background. 

The signal mEs distribution is modeled as the sum 
of two components. The first is a modified Gaussian 
function that, for values less than the mean, has a 



width parameter that scales linearly with the distance 
from the mean. The second component, accounting for 
less than 6% of the distribution, is a threshold func- 
tion [10], which is a phase-space distribution of the form 

"^ES^(1 - ) exp(^(l - ^2"^^ )), with a kinematic 

cut-off at Abeam = 5.289 GeV and one free parameter ^. 
The signal AE distribution is modeled by the sum of 
a Gaussian core with an asymmetric power-low tail [11] 
and a second order polynomial. The parameters of these 
PDFs are determined by fitting to a signal Monte Carlo 
sample. The peak position of the AE distribution is a 
free parameter of the full CP likelihood fit to allow for 
EMC energy scale uncertainties. 

The kinematic variables mEs and AE are corre- 
lated in the B° J/ip K°{Tr°Tr°) and inclusive J/ip 
backgrounds, so two-dimensional PDFs are employed 
for these modes. Variably-binned interpolated two- 
dimensional histograms of these variables are constructed 
from the relevant Monte Carlo samples. 

The mEs PDFs for the BB generic and continuum 
backgrounds are modeled by the threshold function given 
above, and the AE PDFs for these two backgrounds are 
modeled by second order polynomials. The parameters 
for these PDFs are obtained from the BB generic Monte 
Carlo sample and the J/ipf^^^ sample. 

The PDFs used to describe the At distributions of 
the signal and background sources are each a convolu- 
tion of a resolution function TZ and decay time distribu- 
tion V: V{At,aAt) = TZ{St,crAt) ®V{Attrue), where At 
and Attrue are the measured and true decay time differ- 
ences, St = At — Attrue, and a^t is the estimated event- 
by-event error on At. 

For the signal, the resolution function consists of the 
sum of three Gaussian distributions, the parameters of 
which are determined from the Bflav sample, as in the 
B° J/ip measurement [9]. The decay time distribu- 
tion is given by Eq. 1 modified for the effects of B-fiavor 
tagging: 



g-IAtl/r^o 

± 5/ (1 ■ 

TCf (1 



-{(It Awo 



where , (V 



is for a B"(B 



2Wa 
0/- rO\ 



I sm{AmdAt) 
)cos{AmdAt)}, 



(3) 



The 



^j; c ^ ; taggmg mcsou. 

variable Wa is the average probability of incorrectly tag- 
ging a B° as a B° [w^ ) or a B° as a B° [w^ ), and 
Awa = — . Both Wa and Awa are deter- 

mined using the Bflav data sample [3]. We use the values 
Arud = 0.489 ps-i and Tgo = 1.542 ps [12]. 

The PDF used to model the At distribution for the 
B° J/ip li'^^TT^TT^) background, which also includes a 
CP asymmetry, takes the same form as that for signal, 
but with Sj/^x" = sin2/3 = 0.74 [3] and Cj^^x" = 0- 
The parameterizations of the At PDFs for the inclu- 
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TABLE II: Results of the CP likelihood fit, for the full region 
-0.4 < AB < 0.4 GeV and mss > 5.2GeV/c^ Errors are 
statistical only. The global correlation coefficient is 0.f4 for 
Cj/^,rO and O.fS for Sj/^^o. 

Fit results 



C 



S 



Signal AE peak position ( MeV) 

B° J/ij) tt" signal (events) 

B° J/ij) K%{tt° TT°) background (events) 

Inclusive J/j/" background (events) 

BB generic background (events) 

Continuum background (events) 



0.38 
0.05 
-f3.2 
40 
f40 
f09 
52 
97 



0.4f 
0.49 
7.2 

7 

f9 

35 

25 

22 



sive Jlip and BB generic backgrounds each consist of 
prompt and exponential decay components. Decays ap- 
pear to be prompt when particles from the reconstructed 
B are erroneously included in the tagging B vertex. For 
the BB generic background, the prompt and exponen- 
tial components correspond to the cases where the two 
decay products forming the J/^ come from both or just 
one of the B mesons, respectively. The fraction that is 
in the exponential component, the decay lifetime param- 
eter, and the resolution parameters are determined from 
the Monte Carlo simulation. 

The At PDF for the continuum background has only a 
prompt component and the resolution parameter values 
are obtained by fitting the J/tpf^^^ sample. 

The results of the CP asymmetry fit, for all free pa- 
rameters, are shown in Table IT There are 40 ± 7 signal 
events in the total sample of 438 selected events. The 
projection in mEs is shown in Fig. 1. The yields and 
asymmetry as functions of At, overlaid with projections 
of the likelihood fit results, are shown in Fig. 2. Repeat- 
ing the fit with the added constraint Cj/^ j^o = does not 
significantly change the result for Sjf^ j^o . 

The dominant contributions to the systematic errors 
in Cjj^T^o and Sjj^t^o are summarized in Table IIT The 
first class of uncertainties are those obtained by varia- 
tion of the parameters used in the mEs , AE", and At 
PDFs, where the dominant sources are the uncertain- 
ties in the signal AE PDF parameters. A systematic 
error to account for a correlation between the tails of 
the signal mEs and AE distributions is obtained by us- 
ing a two-dimensional PDF. Another contribution stems 
from the impact of EMC energy scale uncertainties on 
the modeling of the 5*^ J/ip K^^n^n^) background. 
An additional systematic uncertainty comes from the 
choice of the binning of the two-dimensional PDFs for 
the B° J/ip K^^TT^TT^) and inclusive J/ip backgrounds. 

In summary, an unbinned extended maximum like- 
lihood fit yields 40 ± 7 signal events and the param- 
eters of time-dependent CP asymmetry for the decay 
B° J/i)TT°: Cji^^o = 0.38 ± 0.41 (stat) ± 0.09 (syst) 
andS'j/^^o =0.05 ±0.49 (stat) ±0.16 (syst). Within the 




Signal + backgrounds 

Sum of backgrounds 



5.24 



5.26 



5.28 



5.3 



tTipg (GeV/c ) 



EIG. 1: Projection in mES for the results of the CP fit, dis- 
played with the added requirement —0.11 < AE < 0.11 GeV. 
In contrast, the CP fit uses the full AE region. In the further 
restricted region niES > 5.27GeV/c , there are 49 data events 
(points), of which about 12 events are fit as background. Here, 
B J/tpKs{7^ TT ) and inclusive J/tp decays contribute to 
the enhancement in the background distribution at large niES- 



Signal + backgrounds ; 
Sum of backgrounds : 




EIG. 2: Distributions of events a) with a B° tag (Ngo), b) 
with a B° tag (A^^o), and c) the raw asymmetry {Ngo — 
Ngo)/{Ngo + Ngo), as functions of At. Candidates in these 
plots are required to satisfy —0.11 < AE < 0.11 GeV and 
niES > 5.27GeV/c . Of the 49 signal and background events 
in this region, 25 have a B tag and 24 have a B tag, with 
fit background contributions of approximately 5 and 7 events, 
respectively. The curves are projections that use the values 
of the other variables in the likelihood to determine the con- 
tributions to the signal and backgrounds. 



Standard Model formulation of CP asymmetries, these 
results demonstrate the possibility, with additional inte- 
grated luminosity, of observing penguin contributions in 
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TABLE III: Summary of systematic uncertainties. 



Source 






Parameter variations 


mES and Ai? parameters 


0.05 


0.13 


Tagging fractions 


0.00 


0.01 


At parameters 


0.03 


0.02 


Additional systematics 


AE-niES correlation in signal 


0.07 


0.08 


EMC energy scale B° J/tP K°s{n°n°) 


0.01 


0.00 


Choice of two-D histogram PDEs 


0.01 


0.03 


Beam spot, boost /vtx., misalignment 


0.01 


0.01 


Total systematic uncertainty 


0.09 


0.16 



B° J/ipTr°. Such a measurement may experimentally 
constrain similar amplitudes in B° J/ipK^. 
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TABLE II: Results of the CP likelihood fit, for the full region 
-0.4 < AB < 0.4 GeV and mss > 5.2GeV/c^ Errors are 
statistical only. The global correlation coefficient is 0.f4 for 
Cj/^,rO and O.fS for Sj/^^o. 

Fit results 



C 



S 



Signal AE peak position ( MeV) 

B° J/ij) tt" signal (events) 

B° J/ij) K%{tt° TT°) background (events) 

Inclusive J/j/" background (events) 

BB generic background (events) 

Continuum background (events) 



0.38 
0.05 
-f3.2 
40 
f40 
f09 
52 
97 



0.4f 
0.49 
7.2 

7 

f9 

35 

25 

22 



sive Jlip and BB generic backgrounds each consist of 
prompt and exponential decay components. Decays ap- 
pear to be prompt when particles from the reconstructed 
B are erroneously included in the tagging B vertex. For 
the BB generic background, the prompt and exponen- 
tial components correspond to the cases where the two 
decay products forming the J/^ come from both or just 
one of the B mesons, respectively. The fraction that is 
in the exponential component, the decay lifetime param- 
eter, and the resolution parameters are determined from 
the Monte Carlo simulation. 

The At PDF for the continuum background has only a 
prompt component and the resolution parameter values 
are obtained by fitting the J/tpf^^^ sample. 

The results of the CP asymmetry fit, for all free pa- 
rameters, are shown in Table IT There are 40 ± 7 signal 
events in the total sample of 438 selected events. The 
projection in mEs is shown in Fig. 1. The yields and 
asymmetry as functions of At, overlaid with projections 
of the likelihood fit results, are shown in Fig. 2. Repeat- 
ing the fit with the added constraint Cj/^ j^o = does not 
significantly change the result for Sjf^ j^o . 

The dominant contributions to the systematic errors 
in Cjj^T^o and Sjj^t^o are summarized in Table IIT The 
first class of uncertainties are those obtained by varia- 
tion of the parameters used in the mEs , AE", and At 
PDFs, where the dominant sources are the uncertain- 
ties in the signal AE PDF parameters. A systematic 
error to account for a correlation between the tails of 
the signal mEs and AE distributions is obtained by us- 
ing a two-dimensional PDF. Another contribution stems 
from the impact of EMC energy scale uncertainties on 
the modeling of the 5*^ J/ip K^^n^n^) background. 
An additional systematic uncertainty comes from the 
choice of the binning of the two-dimensional PDFs for 
the B° J/ip K^^TT^TT^) and inclusive J/ip backgrounds. 

In summary, an unbinned extended maximum like- 
lihood fit yields 40 ± 7 signal events and the param- 
eters of time-dependent CP asymmetry for the decay 
B° J/i)TT°: Cji^^o = 0.38 ± 0.41 (stat) ± 0.09 (syst) 
andS'j/^^o =0.05 ±0.49 (stat) ±0.16 (syst). Within the 
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5.26 
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EIG. 1: Projection in mES for the results of the CP fit, dis- 
played with the added requirement —0.11 < AE < 0.11 GeV. 
In contrast, the CP fit uses the full AE region. In the further 
restricted region niES > 5.27GeV/c , there are 49 data events 
(points), of which about 12 events are fit as background. Here, 
B J/tpKs{7^ TT ) and inclusive J/tp decays contribute to 
the enhancement in the background distribution at large niES- 



Signal + backgrounds ; 
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EIG. 2: Distributions of events a) with a B° tag (Ngo), b) 
with a B° tag (A^^o), and c) the raw asymmetry {Ngo — 
Ngo)/{Ngo + Ngo), as functions of At. Candidates in these 
plots are required to satisfy —0.11 < AE < 0.11 GeV and 
niES > 5.27GeV/c . Of the 49 signal and background events 
in this region, 25 have a B tag and 24 have a B tag, with 
fit background contributions of approximately 5 and 7 events, 
respectively. The curves are projections that use the values 
of the other variables in the likelihood to determine the con- 
tributions to the signal and backgrounds. 



Standard Model formulation of CP asymmetries, these 
results demonstrate the possibility, with additional inte- 
grated luminosity, of observing penguin contributions in 
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